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Abstract. Extensive photoionization model grids for single star Hn regions using a variety of recent state-of-the-art stellar 
atmosphere models have been computed with the main aim of constraining/testing their predicted ionizing spectra against recent 
ISO mid-IR observations of Galactic Hn regions, which probe the ionizing spectra between ~ 24 and 41 eV thanks to Ne, Ar, 
and S fine structure lines. Particular care has been paid to examining in detail the dependences of the nebular properties on the 
numerous nebular parameters (mean ionization parameter U, abundances, dust etc.) which are generally unconstrained for the 
objects considered here. Provided the ionization parameter is fairly constant on average and the atomic data is correct these 
comparisons show the following: 

- Both recent non-LTE codes including line blanketing and stellar winds (WM-Basic and CMFGEN) show a reasonable agreement with 
the observations, although non-negligible differences between their predicted ionizing spectra are found. On the current basis none of the 
models can be preferred over the other. 

- The softening of the ionizing spectra with increasing metallicity predicted by the WM-Basic models is found to be too strong. 

- We confirm earlier indications that the CoStar atmospheres, including an approximate treatment of line blanketing, overpredict somewhat 
the ionizing flux at high energies. 

- Both LTE and non-LTE plane parallel hydrostatic atmosphere codes, predict ionizing spectra which are too soft, especially over the energy 
range between 27.6, 35.0, and 41.1 eV and above. The inclusion of wind effects is crucial for accurate predictions of ionizing fluxes. 

- Interestingly, blackbodies reproduce best the observed excitation diagrams, which indicates that the ionizing spectra of our observed objects 
should have relative ionizing photon flux productions Qe at energies 27.6, 35.0 and 41.1 eV close to that of blackbody spectra. 

These conclusions are found to be fairly robust to effects such as changes of U, nebular and stellar metallicity changes, and 
the inclusion of dust. Uncertainties due to atomic data (especially for Ar) are discussed. We also discuss the possibilities and 
difficulties in estimating absolute stellar temperatures from mid-IR line ratios, or softness parameters defined in analogy with 
the optical 77 indicator n -s=([On] 3726,27 A/[0 m] 4959,5007 A) / ([Sn] 6717,31A/[S in] 9069,9532A). [Arm/n] is found to 
be the only fairly stable T cff indicator. Mid-IR 77's appear to be of limited diagnostic power both empirically and theoretically. 
Finally we have examined which parameters are chiefly responsible for the observed mid-IR excitation sequences. The galactic 
gradient of metallicity changing the shape of the stellar emission is found to be one of the drivers for the excitation sequence of 
Galactic H n regions, the actual contribution of this effect being finally atmosphere model dependent. We find that the dispersion 
of T e fF between different H 11 regions due to statistical sampling of the IMF plus additional scatter in the ionization parameter 
are probably the dominant driver for the observed excitation scatter. 
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galaxies. Because of their intense nucleosynthesis, they pro- 
cess large amounts of material, on very short time scales. 
Furthermore, in addition to type II supernovae, of which they 
are progenitors, massive stars drive the dynamics and energet- 
ics of the ISM through their supersonic massive winds, thus af- 
fecting the subsequent star formation process in their surround- 



1. Introduction 

Despite their paucity, hot massive stars are prominent contrib- 
utors to the chemical and dynamical evolution of their host 
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ing environment. Besides, their strong UV radiative fluxes ion- 
ize the ISM and create H II regions. The ionization structure 
of the latter is therefore, for the most part, controlled by the 
EUV radiation field of their massive stars content. In order to 
determine the properties of H II regions, it is therefore essential 
to understand the physical properties of massive stars and most 
importantly, to constrain their FUV and EUV (H-ionizing con- 
tinuum) flux distribution. Yet, this part of the stellar spectrum 
is generally unaccessible to direct observations and it is cru- 
cial to find indirect tests to constrain it. In this context, nebular 
observations of H II regions combined with extensive grids of 
photoionization models including state-of-the-art model atmo- 
spheres offer the best opportunity to achieve this goal. 

In fact a large number of galactic H II regions have been ob- 
served with the ISO satellite (see e. g lMartfn-Hernandezl2002l 
and references therein). These spectra provide a wealth of spec- 
tral information, through fine-structure lines of ions whose ion- 
ization/excitation threshold are located below 9 1 2 A. The shape 
of the SED in the EUV, and more specifically the number of 
ionizing photons in this region, is directly probed by ratios 
of successive ionization states such as [Arm] 8.98/mi/[Arn] 
6.98/mi, [Nra] 57.3//m/[Nn] 121.8/zm, [Siv] 10.5//m/[Sm] 
18.7/im, and [Nem] 15.5/im/[Nen] 12.8/im. Building line ra- 
tios diagrams for these species that are very sensitive to dif- 
ferent parts of the flux distribution below the Lyman threshold 
allow one to derive informations on the actual spectral energy 
distribution at wavelengths usually unaccessible to direct ob- 
servations. This not only provides valuable informations on the 
physical properties of the H II regions but on their stellar con- 
tent as well. As a matter of fact, it is nowadays often used to es- 
timate the spectral type of the ionizing source of single star H II 
regions, and offers a useful counterparts to more classical tech- 
niques o f typing, based on optical or near-infrared ab sorption 
features (Mathvs 1988; Hanson et al. 1996; Watson & Hanson 
IT99llKaneretall2002l). 

On the other hand, modeling tools to analyze the photo- 
sphere and winds of hot, massive stars with a high level of accu- 
racy and reliability have become available in the recent years. 
In particular, major progress has been achieved to model the 
stellar photosphere and stellar wind in an unified approach in- 
corporating also a treat ment of non-LTE line b lanketing fo r the 
major opacity sources JHil lier ffi: Miller||l998t|Pauldr ach et alJ 
l200lUHubenv & LanJl995ULanz & Huben\l2003albl) . 

The impact of the first generation of atmosphere models 
including stellar winds and non-LTE line blanketing o n neb- 
ular diagnostics was studied by Stasinska & Schaerer (1997) 
using the CoStar atmosphere models of Schaerer & de Koter 
(1997). This study showed already several impro vemen ts with 
respect t o the widely used LTE mo dels of Kuruc zldl99ll). M ore 
recent ly iMartfn-Hernandezl J2002I) : iMartfn-Hernandez et alJ 
(2004)' have investigated the metallicity dependence of the 
spectral energy distribution of O stars and the ionization 
structure of H II reg ions, using the CMFGEN code by 
Iffillier & Miiierl J 19981 . They also compa red the EUV fluxes 
from CMFGEN to thos e of the CoStar jS chaerer & de Koter 
Il997l) and WM-Basic JPauldrach et all200ll) codes. Thev con- 
cluded that different treatment of line-blanketing between 
CoStar on the one hand and WM-Basic and CMFGEN on the 



other hand results in significant differences in the predicted 
EUV SEDs and ionizing fluxes. 

In this context, it is of special interest to investigate how the 
different models available nowadays, compare to each other in 
predicting nebular lines ratios. Similarly, it is of importance to 
test the role that a handful of various nebular parameters might 
have on the line ratios diagrams provided by ISO observations. 
The parameters influencing the ionization structure of a pho- 
toionized region are: 1) the geometry, the density distribution, 
the metallicity of the gas, and the possible absorption of the 
ionizing radiation by dust, 2) any physical quantity affecting 
the shape of the ionizing flux like, for example, the effective 
temperature of the ionizing star, its metallicity, the presence 
of a wind and the characteristics of the latter, 3) the hypothe- 
sis used to model the atmosphere like the number of elements 
taken into account, the treatment of the line-blanketing, etc... in 
summary: the physical ingredients and the related assumptions 
used to model the emitting atmosphere. 

The present paper describes photoionization models per- 
formed with various atmosphere models, separating the effects 
of all these parameters. The paper is structured as follows: The 
various adopted atmosphere models are briefly described and 
compared in Sect. |2] The ionizing spectra from these mod- 
els are then used as input to our photoionization code for the 
calculation of extended grids of nebular models (Sect. [3}. The 
compilation of ISO observations of H n regions is described in 
Sect. 0] In Sect. |3 we compare our photoionization models to 
the observations and discuss the effect of changing parameters 
one by one on the excitation diagnostics. In Sect.|6]we test the 
validity of the different excitation diagnostics and softness ra- 
diation parameters for the determination of T e jf . The discussion 
takes place in Sect.0Our main conclusions are summarized in 
Sect.Ej 



2. O star atmosphere models 

Among the key ingredients for the description of O star model 
atmospheres are the treatment of non-LTE effects, the inclusion 
of stellar winds, and a treatment of line blanke t ing (see e.g. 
[Abbott & Hummer! Il985l iKudritzki et ail 1 19881 iGabler et alJ 
1989). In recent years considerable improvements have been 
made in the modeling of these processes and model grids com- 
puted with various sophisticated atmosphere codes have be- 
come available (see e .g. the recent confer ence on "Stellar at- 
mosphere modeling ". iHubenv et ai1l2003h . For our photoion- 
ization models, we adopt the ionizing spectra predicted from 
five different codes (Kurucz, TLUSTY, CoStar, WM-Basic, 
CMFGEN) briefly described hereafter. With the exception of 
the TLUSTY and Kurucz models, which assume a plane paral- 
lel geometry and thus no wind, all models describe the photo- 
sphere and winds in spherical geometry, in a unified manner. 

Except mentioned otherwise, we have used atmosphere 
models computed for solar abundances: He, C, N, O, Ne, Si, S, 
Ar and Fe being 0.1, 4.7xl0~ 4 , 9.8xl0" 5 , 8.3xl0~ 4 , 5.4xl0~ 5 , 
4xl0~ 5 , 1.6xl0~ 5 , 6.8xl0~ 6 and 4xl0~ 5 resp. 
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Fig. 1. Position in a log(g) versus T e fr diagram of the Supergiant 
models and the Kurucz dwarf models used in this paper. See 
Fig Elfor the line symbols. CoStar models are labeled accord- 
ing to the original naming convention. CMFGEN models have 
the same parameters than WM-Basic models and are not drawn 
here. 

2.1. Kurucz models 

We u se the well-known pla ne parallel LTE line blanketed mod- 
els o flKuruczNl99lUl994l) . Computations were done for mod- 
els with T e ff (and log(g)) between 26 and 50 kK (3.0 and 5.0). 
For stability reasons, the available high T e ff models are re- 
stricted to cases of high gravity. The employed Kurucz models 
are therefore representative of dwarfs rather than supergiants 
mostly considered for the other model atmospheres (cf. below). 

2.2. TLUSTY 

A grid of plane-parallel non-LTE line blankete d mod els based 
on the TLUSTY code of iHubenv & Lanzl i ll 995b has re- 
cently been calculated using a super-level approach and an 
Opacity distribution function or a modified opacity sampling 
jLanz & Hubenvll2003albl) . About 100,000 individual atomic 
levels have been included, for more than 40 ions of H , He , C , 
N , O , Ne , Si , P , S , Fe And Ni , using a superlevel approach. 
For all the models, a standard microturbulent velocity V tul -b=10 
kms~' has been used. The parameter space of the grid covers 
27 500 K < T eff < 55 000 K with 2500 K steps and 3.0 < log(g) 
< 4.75 with 0.25 dex steps. Up to 10 different metallicities, 
from 2 times solar to metal free chemical co mposition, have 
been considered bv lLanz & Hubenvl J2003albh 1 . We extracted 
from this database models with T e g (log(g)) ranging from 30 to 
50 kK (3.0 to 4.0), with solar metallicity. 



ing effects with an opacity sampling method based on Monte- 
Carlo simulations (Schmutz 1991). The impact of these effects 
on the ionizing fluxes and nebul ar diagnostics of H n reg ions 
has been discussed in detail by Stasiriska & Scha ereiHl997l) . 

For our computations we use CoStar models with the low- 
est value for log(g), i.e. models D5, D4, D3, E3, F3, F2 and Fl 
from the CoStar grid of ISchaerer & de Koterl dl997l) . The T eff 
(and log(g)) range from 27 kK (2.9) to 53 kK (4.1). 

2.4. WM-Basic 

The WM-Basic models of P auldrach et alJ J200ll) treat a large 
number of ions in non-LTE and include their line blocking ef- 
fect by means of an opacity sampling technique. The atmo- 
spheric structure is computed from the hydrodynamic equa- 
tions including radiative acceleration from numerous metal- 
lines and cont inua. We used the gr id available on the web 2 and 
described in Pauld rach et alJfeOOll) for Supergiant models with 
Teff (and log(g)) ranging from 30 kK (3.0) to 50 kK (3.9). 

2.5. CMFGEN 

We have constructed spherically symmetric wind mod- 
els, using the non-LT E, comoving frame code CMFGEN 
jHiiTier & Miller 1998). This code solves the radiative trans- 
fer equation, together with the statistical equilibrium equations, 
and line blanketing is self-consistently taken into account, us- 
ing a super-level formulation. The chemical elements included 
in our model calculations are H , He , C , N , O , S , Si and Fe . 
For the 28 ions explicitly treated, a total of 2292 levels dis- 
tributed in 819 superlevels are included, representing 22762 
bound-bound transitions. Atomic data for Feiv and Fe V have 
been slightl y improved, com pared to those first introduced in 
CMFGEN JHillier&Millerlll998l) a nd made consistent with 
those used in TLUSTY. As shown in Bour et et alJ (J2003), this 
was required to get a very good agreement in the determina- 
tion of iron abundances, when fitting lines of these two suc- 
cessive ionization stages in individual O stars in NGC 346, the 
largest H II region in the SMC. The temperature structure is 
calculated from the assumption of radiative equilibrium. The 
atmospheric structure consists of the wind, parameterized by 
the classical /Haw, which is conne cted to hydrostatic laye rs ob- 
tained from the ISA- Wind code of de Kote r et all (11996 ). such 
that at the connecting point both the velocity and velocity gra- 
dient match. We assume a constant Doppl er pro file of 20 km s -1 
for all lines. As shown bv lMartins et al.l J2002I) for dwarfs and 
by additional test calculations this assumption leads to negli- 
gible changes of emergent spectrum. The stellar parameters, 
including the abundances, used to compute the CMFGEN grid 
of supergi ants are identica l to tho se used by WM-Basic and de- 
scribed in Pauldrach^LaDteOJHl)- 



2.3. CoStar 

The CoStar models o flSchaerer & de Koterl ill 9971) include stel- 
lar winds, treat H-He in full non-LTE, and include line blanket- 



2.6. Atmosphere models for Dwarf stars 

The main results of the present paper are obtained for 
Supergiant stars. Nevertheless, we also computed grids of pho- 



1 This grid is available at http://tlusty.gsfc.nasa.gov 



2 http://www.usm.uni-muenchen.de/people/adi/Models/ 
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Fig. 2. Comparison between the 6 stellar atmosphere models: CoStar (solid), WM-Basic (dotted), CMFGEN (dashed), TLUSTY 
(dash dot), Kurucz (dash dot dot) and the Blackbody (long dashes, left panels only), for the same T e ff of 35 kK (upper plots), 
except for CoStar model (see text), and 40 kK (lower plots). The left panels show the Spectral Energy Distribution and the right 
panels show, for any energy E(eV), the number Qe of photons with energy greater than E, relative to the corresponding number 
for the Blackbody emission, all the spectra having the same value for Qo.6- Vertical lines are plotted at 13.6 eV (solid) and 27.6, 
35.0 and 41.1 eV (dotted), corresponding to the ionization potentials of the ions considered in this paper (Ar + , S ++ , and Ne + 
resp.). 



toionization models using Dwarf stellar atmosphere models as 
ionizing spectrum, to check the effect of log(g) on the excita- 
tion of the nebula. In this purpose, WM-Basic D models from 
IPauldrach et alJ ( EoOll) have been used. We have also computed 
CMFGEN models using the same set of parameters than those 
used for the WM-Basic D models. 

The models using Dwarf stellar atmosphere are discussed 
in Sect. 1531 

2. 7. Rebinning of the ionizing spectra 

For subsequent use in our photoionization code NEBU (de- 
scribed in Sect. |5} the different atmosphere models have to 
be rebinned to the energy grid used in this code. The SEDs 
are first converted to the units used in NEBU (number of 
photons/eV/s/cm 2 ). The SED is then interpolated on the NEBU 



grid, such as to preserve the integrated number of photons in 
each energy interval. For most of the energy intervals, the num- 
ber of points in the original stellar atmosphere domain is some 
tens to some hundreds, giving a good accuracy for the rebin- 
ning. Note that despite the low number of points (144) used to 
describe the ionizing spectrum in NEBU, the results are reli- 
able, as the most important quantities are the number of pho- 
tons able to ionize the different ions. The grid points actually 
correspond to the ionization potential of all the ions taken into 
account in the photoionization computation. 

2.8. Comparing the EUV spectra 

Fig. \l\ present all the Supergiant models used in this paper in 
a log(g) versus T e ff diagram. The values for log(g) at a given 
T e ff are very close together, with the exception of the Kurucz 
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models, which have a systematic higher value for log(g), up to 
be even highe r than the value adopted fo r Dwarf models (see 
also Fig. 1 in ISchaerer & de Koterll997l) . 

Fig. |2 illustrates the differences in the SED obtained from 
different atmosphere models after the rebinning procedure de- 
scribed above, for the same T e ff, here 35 and 40 kK, with the 
exception of CoStar model for which no value at 35 kK is avail- 
able in the Supergiant subset of models used here, model D4 at 
32.2 kK is plotted. While the five models agree quite well in 
the domain of energies lower than 20 eV (and very well in the 
optical and IR range, not shown here), their differences can be 
as big as 4 orders of magnitudes just before 4 Rydberg. In this 
paper, we will use IR lines to trace the SED between 27, 35 and 
41 eV (see next Section), where the models differences already 
reach 1 to 2 orders of magnitude. 

Of more interest for the analysis of the behavior of the fine- 
structure lines is the distribution of the ionizing photons at each 
energy. This is quantified by Qe, which is the number of pho- 
tons with energy greater than E, shown in right panels of Fig. 12 
More precisely, the relevant quantity determining the nebular 
structure and properties would be the photon output weighted 
by the photoionization cross section. In the range of energy 
traced by the excitation diagnostics, 27-41 eV, the behavior of 
the four models is very different. We will discuss this further in 
Sect.EUandlO] 



3. Grid of photoionization models 

Extensive grids of photoionization models were computed with 
the NEBU code (|M orisset & Peauignol ll99rjl:IPeauignot et alJ 
l200lUMorisset et all200 2) in order to evaluate in detail the de- 
pendence of the mid-IR atomic fine-structure line emission of 
Galactic Hn regions on the atmosphere models, and the stel- 
lar and nebular properties. Our main aims are a) to derive con- 
straints on the stellar ionizing spectra and b) to examine the ori- 
gin of the observed excitation gradients in (compact) Galactic 
Hii regions. 

In principle nebular emission line properties depend on 
fairly a large number of parameters, namely: 

- the shape of the stellar ionizing spectrum, determined 
(mostly) by the stellar temperature T e ff , gravity, and metal- 
licity Z*; 

- the ionization parameter U(r) = (2i3.6/M>4;rr 2 c. As the ge- 
ometry of the H ii regions modeled here is an empty cavity 
surrounded by an ionized shell, we prefer to use hereafter 
the mean ionization parame ter U = U(r), computed follow- 
ing |Evan^^3o2i^([l985|) at a distance from the ionizing 
star r - r empty + AR/2, where r empty is the size of the empty 
cavity and AR the thickness of the Hii shell 3 . U is essen- 
tially given by the ratio of the ionizing photon density over 
the nebular particle density, i.e. properties of the ionizing 
source and the nebular geometry; 

- the nebular abundances/metallicity Zc, as 

- atomic parameters driving the ionization equilibrium and 
line emissivities, e.g. photoionization cross section, recom- 



bination coefficients (radiative and dielectronic), collisional 
excitation cross sections, etc. 
- other secondary parameters like the presence of dust. 

There is no doubt on the existence of a systematic metallicity 
variation among the observed sources considered below. On the 
other hand, as for most of these objects (compact/ultra-compact 
H n regions) the properties of the ionizing source(s) and their 
geometry are not known, it is imperative to assess the impact of 
all parameters on the observables and to establish that the con- 
clusions drawn from comparisons with observations are robust 
in this respect. The dependence of the emission line properties 
on the above parameters is examined in Sect. [5] with the help 
of photoionization models computed for a wide range of model 
parameters. 

The bulk of "standard" models were computed for the fol- 
lowing cases/assumptions. The ionizing spectra from the five 
atmosphere models described in Sect. 13 and plotted in Pig. El 
plus blackbody spectra are adopted. Stellar T e ff ranging from 
30 to 50 kK were used. This range in T e ff is likely t o describe 
the physical conditions of the sample of H n regions llMorissetl 
2003). For most cases we assume a solar composition for the 
nebular and stellar abundances. Metallicity variations are con- 
sidered in Sect. 15. 51 For each of these stellar atmosphere, series 
of photoionization models were computed for the following 
nebular parameters. We set the electron density to 10 3 cm -3 , 
one order of magnitude belo w the lowest critical density of th e 
lines under consideration (cf . iMartfn-Hernandez et al1Eo 02a). 
An empty cavity of radius 3xl0 17 cm is assumed. The luminos- 
ity of the ionizing star is adjusted to lead to a constant num- 
ber of Lyman continuum photons (Q13.6 = 1.5 10 49 s -1 ) corre- 
sponding to an ionization parameter log(£7) = -1.5. Additional 
models quantifying the effect of variations of U are presented 
in Sects.ESEl andIO 

The effect of dust can be included in the photoionization 
computation, with two different optical properties correspond- 
ing to graphite or astronomical silicate (see Sect. 15. 6> . 

The observables predicted from these extensive model grids 
will be compared to observations in Sect. 15. II 

4. The ISO H 11 regions catalogs 

Infra-red spectra between 2.3 and 196 fim were taken from a 
sample of 43 compact H 11 regions using the two spec trometers 
(SWS and LWS) on board ISO JPeeters et alJl2002l) . Details 
about the data reduction and a first ana lysis of the ionic lines in 
terms o f abundances can be found in Martin-Herna ndez et alJ 
(2002a). Error bars on the lines intensities are within 10% to 
20%. Note that a detailed study of one source has been achieved 



in Morisset et al. ( 


2002 


. 




Giveon et al. ( 


2002b) published a catalog of 112 Hii re 



3 For U derived at the Strbmgren radius without empty cavity one 
has U oc (QN e e 2 ) lli , with the filling factor e. 



gions observed by ISO SWS spectrometer. Some of the sources 
are common with the Martin-Hernandez et alJ J2002al) catalog. 
The two catalogs are very coherent in t erms of line intensities, 
as concluded by iGiveon et alJ d2002ah . and are therefore in- 
cluded in our analysis. The effect of local and interstellar atten- 
uation, even if lower in the IR range used in this work than for 
the optical domain, can be important and need to be corrected 
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for. We correct the observed line intensities fro m the redden- 
ing usin g the extinction law described in Tab. 2 o flGiveon et alJ 
d2002d) . 

In the SWS and LWS spectral domain, 4 fine-structure line 
ratios are sensitive to the ionizing flux distribution: [Arm] 
8.98/mi/[Arii] 6.98yum, [Nra] 57.3/im/[Nn] 121.8/mi, [Siv] 
10.5//m/[Sin] 18.7yum, and [Nem] 15.5/zm/[Nen] 12.8jum, 
hereafter [Ariu/ii], [Nm/ii], [S iv/iii], and [Nem/ii] respectively. 

The excitation ratio implying Nitrogen lines will not be 
used in the next discussion, since: 1) the two Nitrogen lines 
are observed by LWS spectrometer, with a larger aperture 
size than the SWS: some Nitrogen emiss ion can ar i se from 
regions not seen in the other lines; 2) iGiveon et alJ (12002b) 
have observations only with SWS and then the number of 
observational constraints strong ly decrease when using only 
Martin-Hernandez et al. (2002 a) data; 3) the critical densities 
of the Nitroge n lines are very low compared to the one of the 
other lines (see Martin-Hernandez et al. 2002a) and will not be 
emitted by medium density gas which can still emit the other 
lines, and 4) the ionization potential of N ++ is very close to 
the one of Ar ++ (29.7 and 27.6 eV resp.), so the main conclu- 
sions regarding the 30 eV energy domain will be obtained from 
Argon lines. 

Depending on the element, the number of sources for which 
we have finite value for the corrected excitation ratios is rang- 
ing from 45 to 5 1 . Error bars on the line intensities are approx- 
imately 10 to 20%. 

5. Excitation diagnostics 

Figs.[3]to[5]show the main results of the photoionization models 
using different atmosphere models for the ionizing star (lines), 
and t he deredenned observed values, for the two merged c ata- 
logs jMartm-Hernandez et al.ll2002at IGiveon et alJl2002bl) . As 
we consider 3 diagnostic ratios, 3 plots can be drawn. The mod- 
els obtained with T e ff = 35 and 40 kK are symbolized by filled 
diamonds and open squares respectively. The open diamond 
indicates a CoStar model at 32.2 kK as no model at 35 kK is 
available. 

In principle the position of a model in such diagrams de- 
pend on all the following parameters: the hardness of the stellar 
SED (parametrized here for each set of model atmospheres by 
T e ff for a fixed stellar metallicity) and the main nebular param- 
eters, i.e. the ionization parameter U and the nebular composi- 
tion. Let us consider first the case of constant (solar) metallic- 
ity and constant V (but see Sects. l5~2l 15. 4l and l5.5> . In this case 
the location of a point on such diagnostic diagrams depends 
only on a) the global excitation of the gas and b) the "slope" of 
the ionizing photon distribution between the two corresponding 
ionization potentials. For constant U and a given set of atmo- 
sphere models the excitation (a) is determined by T e g . In other 
words, when T e ff increases, the number of ionizing photons at 
all the energies traced by the observed ions increases, and the 
points in the excitation diagrams will essentially move along 
the diagonal (y = x) direction. Note, however, that different 
atmosphere models with the same T e fr predict fairly different 
absolute positions in these plots. This simply reflects the dif- 
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[SIV]/[SIII] (35.0eV) 

Fig. 3. Deredenned observed values for the excitation sensitive 
line ratio [Ar m/ii] versus [S iv/m] (the corresponding ionization 
potentials are also given). Source with a galactocentric distance 
lower than 7 kpc are symbolized with a +, otherwise with an 
X. Results from the photoionization model grid are line plotted 
using the same codes as in Fig.|2] The plot have been done such 
as the lowest ionization potential (indicated in braces) is always 
on the y-axis. Models obtained with 35 and 40 kK stars are 
shown using filled diamonds and empty squares respectively 
(except for CoStar model at 32.2 kK, empty diamond, see text). 
The y=x line is also drawn. 
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[Nelll]/[Nell] (41.1eV) 

Fig. 4. Same as Fig. [5] for the excitation sensitive line ratio 
[Arm/ii] versus [Neiu/n] 

ferences in the predicted number of ionizing photons above the 
relevant energy (cf. Fig. [3}. 

For a given line ratio the other line ratios depend to first or- 
der on the "slope" of the ionizing spectra (b). More precisely, 
the relevant quantity is the slope of the cumulative number of 
ionizing photon flux Qe between the corresponding ionization 
potentials (see right panels in Fig. |2}. For example, TLUSTY 
and Kurucz models show in general the softest spectra (i.e. 
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Fig. 5. Same as Fig. [3] for the excitation sensitive line ratio 
[S iv/m] versus [Ne m/ii] 



steepest slopes) between 27.6 and 41.1 eV. For a given [Ne m/ii] 
these models therefore show the highest [Arm/ii] values. 

For the assumptions made here (constant U and metallicity) 
every location of the model results in the three excitation dia- 
grams can be approximately understood in terms of the ioniz- 
ing photon distributions Qe. The correspondence is not always 
exact, as some competitive processes take place in the use of 
the ionizing photons, but the overall trends can be simply un- 
derstood from the shapes of the spectra. Other additional as- 
sumptions (e.g. on the luminosity class of the exciting sources, 
the presence of dust, and uncertainties of the atomic data) also 
affect the predicted excitation diagrams. These effects are dis- 
cussed below. 

The observed excitations, correlated between the three ex- 
citation ratios [Arm/ii], [S rv/m], and [Ne m/n], can be decom- 
posed into two components: an excitation sequence showing 
a global increase of the excitation ratios over ~ 2 orders of 
magnitude, following to first order a trend parallel to y = x in 
the excitation diagrams, and a superposed excitation scatter 
of typically ~ 0.5 - 1 order of magnitude around the mean ex- 
citation (cf. Figs. [3] to [3 but see also Figs. 1191 and 1201 where 
excitations versus metallicity are plotted). 

5.1. First comparison with the observations 

The zero-th order trend of the observations plotted in Figs. [3] 
to[5]is reproduced by the models: the excitation of the ionized 
gas, traced by the X' +1 /X ! ratios, are well correlated. A T e ff 
sequence from ~ 30 to 45 kK succeeds in reproducing the entire 
range of gas excitations. Note however, that, as discussed below 
(Sect, I7.2L this does not imply that the ionizing stars of our 
objects indeed cover this range of T e ff . 

Fairly large differences are found in the predicted excita- 
tion diagnostic diagrams (Figs. [3] to [5} when using different 
atmosphere models. As expected from the intrinsic SEDs, the 
largest differences are found in Fig.0] which traces the largest 



energy domain (~ 28 to 41 eV) corresponding to the [Arm/ii] 
and [Nem/ii] ratios. 

When taken literally (i.e. assuming a fixed constant value of 
U for all atmosphere model sets and a fixed solar metallicity) 
Fig s -|3jto[5]indicate m e following concerning the shape of the 
ionizing spectra. 

1 . Both recent codes, WM-Basic and CMFGEN, show a rea- 
sonable agreement with the observations. Given their differ- 
ent behavior in the three excitation diagnostics depending 
on luminosity class (cf. Sect. 15. 3> . and other remaining un- 
certainties discussed subsequently, it seems quite clear that 
none of the models can be preferred on this basis. 
Despite these similarities we note, however, that an impor- 
tant offset is found in the excitation ratios predicted by these 
codes for a given absolute value of T e jj (cf. Sect. 16. 1> . 

2. Both plane parallel hydrostatic codes (Kurucz, TLUSTY) 
predict spectra which are too soft, especially over the en- 
ergy range between 27.6 and 41.1 eV and above. For 
the Kurucz LTE models this problem is just a manifes- 
tation of the well-known "Ne m" problem jRubin et alJ 
1988; ISimnson etalJ 1 199.4 IStasiriska & Schaererl ll 9971) 
This problem is persistently found with the non-LTE 
TLUSTY models. Although not completely clear, this "soft- 
ness" is likely due the neglect of win d effects which are 
known to alter the ionizing s pectrum ("cf. lGabler et alJ 1989; 
Stasi hska & Schaererl 1 1997|) albeit in fairly complex way 
involving line blanketing from large numbers of metal- 
lines. 

3. As already found in other investigations (cf. lOev et all 
l2000l ISchaerejl2000l and references therein) we see that 
the CoStar models (showing the hardest spectra among the 
ones considered here) overpredict somewhat the ionizing 
flux at high energies. This likely overestimate of the ioniz- 
ing flux at high energy (cf. below) must be due to the ap- 
proxi mate and incomplete treatment of line b lanketing (see 
also lCrowther et alll999tlMartm-Hernandez et al.l2004 . 

4. Interestingly blackbodies reproduce best the observed ex- 
citation diagrams, which indicates a posteriori that the ion- 
izing spectra should have relative ionizing photon flux pro- 
ductions (Qe at energies 27.6, 35.0 and 41.1 eV) close to 
that of blackbody spectra. This will be discussed in more 
detail in Sect.|7] 

We note that the results concerning CMFGEN, WM-Basic, 
CoStar, and Kurucz m odels presented here c onfirm and sup- 
port those presented bv lStasihska et alJ J2002I) in terms of po- 
sition of the models between each other in the [ Ar m/n] versus 
[Ne m/n] excitation diagram and distance to the H 11 regions ob- 
servations (cf. below). 

5.2. Main dependences: stellar temperature, ionization 
parameter, metallicity 

For clarity it is useful to discuss first the dependences of the 
excitation diagnostics on the main parameters, i.e. the stellar 
temperature, ionization parameter, and metallicity. This is illus- 
trated here somewhat schematically for the case of the [Ne m/ii] 
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versus [Arm/ii] diagnostic. Qualitatively the same results are 
found for the other excitation diagrams. 

An increase of the stellar temperature T e ff or ionization pa- 
rameter U or a decrease of the metallicity all lead overall to 
a higher excitation of the nebula, which e.g. manifests itself 
by larger [Nem/n] and [Arm/n] line ratios. However, although 
both line ratios change in similar ways, their effect is distin- 
guishable to some extent. This is illustrated in Fig. |6] which 
shows for blackbody (and WM-Basic, see Sect. 15. 5> spectra the 
implied shift in the [Neiii/n] versus [Arm/ii] excitation diag- 
nostics due to a change of T e ff , U and Z eas (consistent changes 
of both the stellar and nebular metallicity are discussed in 
Sect. 15.51 . From this figure we see that T e ft and U variations 
are not completely degenerate (i.e. "parallel"). 

Also, an increase of the nebular metallicity leads to a (very) 
small decrease of the excitation diagnostics, quite parallel to 
the variations induced by changing U. The effect of changing 
coherently the metallicity of both the H n region and the star 
(which is principally acting on excitation diagnostics) are con- 
sidered in Sect. 15.51 

The effect of continuum absorption by dust inside the H n 
region on the excitation diagnostics is also parallel to changes 
of T e ff, and is discussed in Sect. 15.61 Although quantitatively 
these variations depend e.g. on the adopted SED (and of the 
point in the T e ff-t/-Z space chosen to compute the partial 
derivatives traced by the arrows in Fig.[6j, these qualitative dis- 
tinctions remain valid for the entire parameter space considered 
in the present paper and will be useful for the discussions be- 
low. 

Note that earlier investigations have co nsidered that 
change s of U are completely dege n erate wi th T e ff JGiveon et alJ 
l2002bl iMartfn-Hernandez et alJ l2002bl) 4 . or have sim- 
ply assumed an arbitrarily fixed, constant value of U 
(Mar tin-Hernandez et all2004 . 

As apparent from Fig. |6] [Ariii/n] varies little with U in 
comparison with other mid-IR excitation ratios. This behavior 
is easily understood, as explained in the following brief digres- 
sion. 

In T e ff and U domain used in this work, the [Arm/ii] ra- 
tio is controlled by the Helium equilibrium: the IPs of Ar + and 
He are closed together (27.7 and 24.6 eV resp.); in this energy 
domain the photon dominant predator is He . The Ar ++ region 
(Ar + ) is then cospatial with the He + (He ) region (Some Ar + 
can also be present in the He + region, depending on the ioniza- 
tion parameter). The H n regions modeled here are all radiation 
bounded, the size (and the emission) of the He + and Ar ++ re- 
gion is mainly proportional to Q24.6, while the size of the He 
and Ar + region is controlled by the size of the H n region re- 
moving the He + region. [Arm/11] is then mainly controlled by 
Q24.6/Qi3.6- The previous argumentation is valid only if the re- 
combination of Ar ++ remains quite small, which is not the case 
when strong dielectronic recombination occurs. In this case, the 
Ar + region penetrates inside the He + region, and the [ Ar m/11] 
is decreased (see Sect. 15.71 for the effects of dielectronic re- 



For their grids the ionization parameter U defined at the 
Stromgren radius vary actually by ~ 0. 1 (supergiant grid) and 0.3 dex 
(dwarfs). 
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Fig. 6. Increase of the excitation diagnostics [Arm/n] versus 
[Nem/n], when log({7) is increased from -2.6 to -0.8 (solid ar- 
rows), when the metallicity of the gas is increased from half 
solar to twice solar (dashed arrows), and when T e ff is increas- 
ing from 35 to 40 kK (dotted arrows). Upper set of arrows are 
for BlackBody models, and lower set for WM-Basic models. In 
the latter case the metallicity is changes coherently for both the 
gas and the star. Codes are the same than in Fips. l2|3l 
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Fig. 7. Correlation between Hei 5876A/H/? and [Arm/n]. The 
line styles are the same as in Fig. |3 Arrows have the same 
meaning as in Fig.|6j for WM-Basic models only. 



combinations of Ar ++ .) Nevertheless, using atmosphere mod- 
els instead of BlackBody leads to a more important increase of 
[Arm/n] while increasing U, as seen in Fig. [6] 

The extreme correlation between He + /H and Ar ++ /Ar + can 
be verified in Fig. where Hei 5876A/H / 3 versus [Arm/n] 
is plotted, for all the atmosphere models. While the Hei 
5876A/Hg ratio saturate at a value between 0.1 and 0.2, the 
[Arm/n] excitation diagnostic still evolve with T e ff. This will 
be discussed further in Sect. 16. II 
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Fig. 8. Comparison between energy distribution (same as right 
panels of Fig.|2) between Supergiants (light curves) and Dwarfs 
(bold curves) of WM-Basic (dotted) and CMFGEN (dashed) 
models, all at 35 kK. 



100.00 



1 0.00 = 



1 .00 r 



0.01 



p I I 


1 1 1 = 




Lai _ 






.... 1/ 1 1 





0.01 0.10 1.00 

[Nelll]/[Nell] (41 ,1eV) 



10.00 



Fig. 9. Comparison between excitation diagnostic [Ar m/n] ob- 
tained with Supergiants (light curves) and Dwarfs stars (bold 
curves), for WM-Basic (dotted) and CMFGEN (dashed) atmo- 
sphere models. Models at 35 and 40 kK are shown by filled 
diamonds and empty squares respectively. 

5.3. Comparing results of Supergiant and Dwarf stars 

The results shown in previous Sections are obtained for 
Supergiant stars. Models were also performed using Dwarfs 
stars (see description in Sect. 12. 6> . The decrease of log(g) 
changes the shape of the ionizing radiation, as shown in Fig.[H] 
where Q E is shown for Supergiants and Dwarfs atmosphere 
models obtained using CMFGEN and WM-Basic, all at 35 kK. 
The main effect of increasing log(g), observed on the two mod- 
els, is to decrease strongly Qe at 41 eV, and to increase the 
slope of Qe between 27 and 41 eV. This overall hardening of 
the ionizing flux for stars with lower gravity is due mostly to an 
increased ionization in the continuum forming layers, the lat- 



ter effect resulting from the increased wind density (mass-loss 
rate). 

Fig. [9] shows the excitation diagram [Nera/n] versus 
[Arm/ii] performed using the models described above, com- 
paring the Supergiant (light curves) and Dwarf (bold curves) 
results for both CMFGEN (dashed) and WM-Basic (dotted). 

As expected from the increased hardness of the ionizing 
fluxes for supergiants, the use of dwarf atmospheres leads 
in general to an excitation decrease which is more important 
at the highest energies (i.e. [Nem/n] decreases more rapidly 
than [ Ar m/ii] ) 5 . Overall the differences between supergiant and 
dwarf spectra do not importantly affect our conclusions. 

5.4. Effect of changing the mean ionization parameter 

U 

Given that the ionizing sources and the nebular geometry of 
the observed objects are essentially unknown and U therefore 
as well, it is important to examine how robust the above results 
are with respect to changes of U . E.g. is it possible to reconcile 
the discrepant predictions using the Kurucz and TLUSTY at- 
mosphere models (i.e. to increase the predicted [Nem/ii] ratio 
of a given [Ann/11]) by invoking a larger ionization parameter 
toward the low excitation end of the observed sequence? As 
shown in Sect. 15.21 and by detailed model calculations, varia- 
tions of U imply changes nearly parallel to the "standard" se- 
quences for the Kurucz and TLUSTY atmosphere models, sim- 
ilar to the case of WM-Basic models shown in Fig.|6] 

We can therefore quite safely conclude that variations of 
U cannot reconcile the Kurucz and TLUSTY atmosphere mod- 
els with the observation. However, for the other model atmo- 
spheres showing more "curved" predictions in the excitation 
diagrams, variations of U might be invoked to improve/alter 
the predicted sequences for constant U . 

5.5. Effect of adopting a metallicity gradient, for both 
the stellar atmosphere and the ionized gas 

It is well known that the metallicity decreases i n the Galaxy 
when t he distance to the center incre ases (e.g. iGiveon et alJ 
2002at lMartfn-Hernandez et al.l2002al and references therein). 
The metallicity varies approximately by a factor 4 from the cen- 
ter out to 15 kpc, where the most external regions used in this 
work are located (cf. Fig. 1191. 

Fig-EJshows the effect of changing the metallicity Z coher- 
ently in the ionizing star's SED computation and of the ion- 
ized gas in the photoionization computation, from half solar 
to twice solar. A metallicity increase in the atmospheres leads 
to a stronger blanketing with t he effect of softening the EUV 
spectra of early type stars (e.g. lPauldrach et al]l200ll) . leading 
thus to a lower nebular excitation. As seen by comparing the Z- 
arrows in Fig.[6]and by additional test calculations, the increase 
of the nebular abundance plays a minor role in the resulting ex- 
citation shift. 



5 The WM-Basic model at 35 kK shown here is an exception to this 
trend, however. 
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In fact the WM-Basic models used here show that the ioniz- 
ing spectra soften too strongly with increasing metallicity, lead- 
ing to a stronger reduction of [Ne m/ii] compared to [ Ar m/ii] . 
This results in a progressive shift away from the observed se- 
quence t oward higher metallicity. This discrepant trend, also 
found bv lGiveon et alJ J2002bl) . was actually used by these au- 
thors to argue that the observed excitation sequence was mostly 
driven by T e ff variations. However, as abundances of these 
sources are known to vary by approximately the same factor 
as the Z variations considered here for the WM-Basic models, 
metallicity cannot be neglected. Therefore the discrepancy be- 
tween the observations and the expected changes of [Ne m/n] 
and [Arm/n] show that the predicted softening of the WM- 
Basic ionizing spectra with metallicity at high energies (^41 
eV) is probably incorrect. Alternate solutions to this puzzle in- 
clude postulating a increase of U toward higher Z (cf. above), 
or processes currently not accounted for in the WM-Basic mod- 
els altering the high energy part of the SED (cf. Sect. 15.81 . or 
changes in atomic physics parameters (cf. Sect l5.7> . 
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5.6. Effect of the dust 

The effect of the presence of dust inside an H n region is firstly 
to decrease the global amount of ionizing photons from the 
point of view of the ionized gas, the effect being then to reduce 
the ionization parameter. On the other hand, the efficiency of 
dust in absorbing of the ionizing photons inside the H n region 
decreases with the energy of the photons after about 18 eV (e.g. 
lMathislll985l lAannestadlll989l) . the global effect being to in- 
crease the excitation of the gas when increasing the amount of 
dust, for a given i onization parameter. As already pointed out 
by iMorisset et al.l J2002I) for the case of ISO observations of 
G29.96-0.02, if dust is present in H n regions, quite the same 
excitation of the gas will be recovered using an higher ion- 
ization parameter and a lower T e ff. As illustration, inclusion 
of Graphite and Astronomical Silicate dusts, in proportion of 
2.5xl(T 3 relative to hydrogen, for each type of dust, leads to 
an increase of all the excitation diagnostic ratios by a factor 
close to 2. The excitation increase due to dust is found to be 
"parallel" to a T e ff increase to reasonable accuracy, when keep- 
ing U constant. 



5.7. Effect of the atomic data 

Could uncertainties in the atomic data affect the results ? 
Indeed, from Fig. [3] to |5] we could suspect the Argon ioniza- 
tion equilibrium to be wrong, favoring the emission of [Arm]. 
This could e.g. be due to an overestimation of the ionizing flux 
at 27.6 eV with respect to higher energies, or to a systematic er- 
ror in the observed intensities of one of the two lines involved in 
the [Arm/ii] ratio ([Arm] 8.98/im is affected by silicate band). 
However, we cannot exclude also the effect of atomic data in 
the photoionization computation. Collisional rates are gener- 
ally believed to be accurate within 20%, while our knowledge 
of recombination coefficients are less probant. Dielectronic re- 
combination coefficients for the elements of the third row of the 
periodic Table are poorly known, and even the new computa- 
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Fig. 10. Variation of the excitation diagnostics [Ne m/n] versus 
[Arm/n] for the same atmosphere models (here WM-Basic and 
CMFGEN, dotted and dashed thin lines respectively), multiply- 
ing the effective recombination coefficient for only Ar ++ by 10 
(bold lines). 



tions done today are usually only for the first and second rows, 
corresponding to highly charged ele ments of the third rows , 
which is not the case for Ar ++ (see e.g. Savin &Lamingl2 002). 
Dielectronic recombination coefficients have been comp uted 
for les s charged third row elements by e.g. Maz zotta et alJ 
but only for high electron temperature (coronal gas), 
which is not the case for H n regions. Very recently, new di- 
electronic recombination rates have been computed for Ar ++ 
(Zatsarinny et al., 2003, private communication), but the re- 
sults these authors obtain have still to be checked (dielectronic 
recombination rates reach values as 10 3 times the classical re- 
combination rates for electron temperature closed to 10 4 K!). 

To simulate the effect of dielectronic recombination and 
charge transfer reactions we have multiplied the classical re- 
combination coefficient for Ar ++ by factors up to 20. Fig. [TO] 
shows the effects of multiplying this coefficient arbitrarily by 
10, on the excitation diagnostics [Arm/n] versus [Nem/ii]. The 
figure shows results using WM-Basic and CMFGEN models, 
but the same effect can be observed with any of the atmosphere 
models used in this paper. 

The increase of the recombination coefficient improves 
overall the agreement with the observed sequence, although 
new discrepancies appear now at the high excitation end. 
However, no dependence of the dielectronic recombination co- 
efficient on the electron temperature T e has been taken into ac- 
count here. As T e is known to vary along the excitation se- 
quence this could in fact "twist" the global shape of the pre- 
dicted excitation sequence. Currently both the exact "direction" 
and importance of this effect remain, however, unknown. Note, 
that the uncertainti es due to atomic d ata of Ar were already 
pointed out by IStasihska et alJ J2002I) . We join these authors 
in encouraging atomic physicists to improve our knowledge of 
such data. 
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Fig. 11. Variation of the excitation diagnostic [Ar m/ii] accord- 
ing to the T e ff of the ionizing star. Codes are the same as in 
Figs. 12131 The set of arrows have the same meaning as in Fig. 
|6j for WM-Basic models only. 

5.8. Effects not included in atmosphere models 

What could be the limitations of the most sophisticated at- 
mosphere models currently available, capable of altering the 
excitation diagnostics discussed here? Although an exhaus- 
tive discussion is obviously not possible, one can suspect one 
major process, namely the presence of X-rays, to alter in a 
non-negligible way the ionizing spectra of O stars. This has 
been shown clearly by Macfarla ne et alJ dl994l). an d has been 
discussed later e.g. bv ISchaerer & de Koterl dl997h . The rela- 
tive importance of X-rays compared to normal photospheric 
emission is expected to increase for stars with weaker winds 
and toward later spectral types. In late O types their contri- 
buti on can be non-ne gligible down to energies sS 30 eV, see 
e.g. iMacfarlane et alJ dl994 and a model at T eff =30 kK by 
Pauldrac h et alJ d200ll) . with obvious consequences on nebular 
diagnostics. Regrettably few models treating the X-ray emis- 
sion in O stars exist, their impact on the overall emergent spec- 
trum including the EUV has hardly been studied with complete 
non-LTE codes including winds and blanketing, and their de- 
pendence with stellar parameters (wind density, stellar temper- 
ature, even metallicity?!) remains basically unknown. 

For now, we can only qualitatively expect the inclusion of 
X-rays to harden the ionizing spectra, probably down to the en- 
ergy range probed by (some) mid-IR diagnostics. While this 
could in principle improve some difficulties observed by the 
CMFGEN and WM-Basic models (e.g. increasing [Nem/ii]) 
their precise effect remains open. 

6. T e ff diagnostics and "second-order" diagnostic 
diagrams 

6.1. Using excitation diagnostic to determine T ef f 

Since the excitation of the gas increases with T e ff , it is tempting 
to infer stellar temperatures from excitation diagnostic ratios. 
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Fig. 12. Variation of the optical excitation diagnostic He i 
5876A/H/J with T e ff for supergiant stars. Codes are the same 
as in Figs. 12131 The set of arrows have the same meaning as in 
Fig-El f° r WM-Basic models only. 
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Fig. 13. Variation of the optical excitation diagnostic [O m] 
5007A/[On] 3727,29 A with T e ff for supergiant stars. Codes are 
the same as in Figs. 12131 The set of arrows have the same mean- 
ing as in Fig. [6] for WM-Basic models only. 

However, such an approach is intrinsically highly uncertain, 
as the nebular excitation is also strongly dependent on other 
parameters (see Sect. such as the ionization parameter U , 
which remain in mo st cases poorly known , cf. pvlathisl dl982l) 
for optical lines and ISchaerer & Stasihska d 19991) for mid-IR 
ratios. These cautionary remarks should be kept in mind when 
e.g. using single line ratios or even several lin e ratios (e.g. 
i Takahashi et al ll2000UOkamoto et al 1200 1 U2003h . but see also 
Morisset (2003), to estimate stellar properties of individual ob- 
jects from nebular observations. Tailored photoionization mod- 
els including numerous constraints can lead to substantially dif- 
fere nt results and sh ould clearly be the preferred method (see 
e.g. Moriss etet al.l20 02). 



12 



C. Morisset et al.: EUV fluxes of massive stars 



For illustration we show in Fig.[^the dependence of the 
[ Ar m/n] excitation ratio on T e g for a fixed U and metallicity. 
Other mid-IR excitation diagnostics show similar behaviors as 
can already be seen from various figures above, except that 
their dependence upon U are higher than for [Arm/n], as al- 
ready discussed in Sect. 15.21 The discussion of T e g and U de- 
terminations using mid-IR excitation diagno stics an d the Hn 
regions metallicities is developed in Morissetl (12003). 

The most important conclusion from Fig.^Jis the impor- 
tant difference of the excitation of the gas ionized by WM-Basic 
(dotted line) and CMFGEN (dashed line) stars, for the same T e g 
and U, even if the two types of models are showing the same 
behavior in Fig. [3] to |3] In other words, taking for example a 
value of 10. for [Arin/n], we can determine a T e ff of 39 kK 
using CMFGEN and a value of 45 kK using WM-Basic. This 
behavior can easily be understood when comparing the Qe dis- 
tribution, as shown in Fig.[2]and discussed in Sect. 12. 81 

Two of the classical ways to constrain T e ff from opti- 
cal observations are to use the Hei 5876A/Hg ratio (e.g. 
iKennicutt et ai1l2000l) or t he [Oiii] 5007A/[On] 3727,29Ara- 
tio JPors & Copettill2003l) . The predictions for He i 58761/^ 
are shown in Fig.^] As for [Arm/ii], this line ratio is fairly 
independent of U. This diagnostic line ratio saturates above 
T e ff <: 40 kK, when Helium is completely ionized to He + . Note 
that even among WM-Basic and CMFGEN, which treat very 
similar physics, some differences in this T e ff indicator remain. 
Furthermore note that the predicted Q24.6/Q13.6 and hence Hei 
5876A/H fl vary non-negligibly between dwarfs an d supergiants 
(see e.g. Paul drach et al.ll2001tlSmith et alJl2002l) . The predic- 
tions for [O m] 5007A/[On] 3727,29 A are presented in Fig.fBl 
As for the He 1 5876A/H J g ratio, the results differ strongly from 
one atmosphere model to another one. The Oxygen excitation 
diagnostic is also strongly sensitive to U and Z, the effect of 
Z being slightly more important than found by Dors & Cooetti 
il2003h . but these authors used WM-Basic dwarfs atmosphere 
models. 

We can conclude from the above examples that any de- 
termination of T e ff based on diagnostic ratios can be reliable 
only if the metallicity of the star is coherently taken into ac- 
c ount, and if U is also determined at the same time, as shown 
in lMorissetlll2003l) . 

6.2. Radiation softness parameters r\ 

Following lVilchez & Paeell dl988t) . radiation softness parame- 
ters can be defined combining the excitation diagnostic ratios, 

[Arm/11] 



namely: j] Ar - Ne 
ratios. 



[Neni/ll] 



To first order (but see the discussion on [Ann/11] in 
Sect. 15. 2> an excitation ratio X' +1 /X' depends on the ionization 
parameter U and the hardness of the ionizing radiation and is 
given by 



X 

~x 



oc U - 



u 



QE(X j ] 



C m r v dv Q 



(1) 
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Fig. 14. Variation of the radiation softness parameter T]Ar-Ne 
versus [Arm/11]. The codes are the same than in Figs. 12131 Note 
that the model predictions for both axis depend also on the as- 
sumed ionization parameter U. Arrows as in Fig.il for WM- 
Basic models. 
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Fig. 15. Variation of the radiation softness parameter r\Ar-Ne ac- 
cording to the T e ff of the ionizing star. Codes are the same as 
in Figs. 12131 Arrows have the same meaning as in Fig. |6] for 
WM-Basic models. 



and so on for the other diagnostic r/x-y 00 



Therefore one has: 

Qe(X') 



(2) 



where E(X') is the ionization potential of the ion X'. r\ 
is thus in principle independent of the ionization parame- 
ter U, and a measure of the "slope" of the ionizing spec- 
trum between the ionization energies E(X') and E(Y') re- 
spectively 6 . Therefore such quantities - used so far for opti- 
cal lines only - have often been thought to be good estima- 



13.6 



6 However, note that e.g. any rj involving [Arm/11] will likely de- 
pend on 0, as [Arm/11] itself depends already little on U 
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Fig. 16. Variat ion of the radiation softness parameter rjs-o (de- 
fined by Vilchez & Pagel ( 1988), see note[8}, according to the 
T e ff of the ionizing star. Codes are the same as in Figs. 12131 
Arrows have the same meaning as in Fig. |6] for WM-Basic 
models, for WM-Basic models only. 

tors of stel lar effective temperatu res, see e.g. |GarnetJri 989): 
IKennicutt etaPfcOOOll but also lSkillmanl dl989l) : lBresolin et alJ 
( 1999J; lOev et al.l (I2002h . For this reason we here explore 
whether the observed mid-IR fine structure lines offer, both 
from the empirical and theoretical standpoint, such a diagnostic 
power. 

Fig. 1141 shows the variation of rj^e-Ar versus [Arm/n]. No 
correlation between these two observables is found. While a 
galactic gradient is found fo r [Ar m/n] JMartm-Hernandez et alJ 
2002a; Giveon et al. 2002b), r\Ne-Ar and none of the other mid- 
IR 77's which can be constructed show a galactic gradient. This 
can also be seen from the correl ations between the vario us ex- 
citation ratios plotted by Martin-Hernandez et al. ( 2002a), their 
Fig. 10. Already this finding indicates empirically that mid-IR 
softness parameters do not carry important information on the 
stellar ionizing sources. 

How do the models compare with the observed softness pa- 
rameters? Given the considerable spread between photoioniza- 
tion models using different stellar atmospheres and the vari- 
ous discrepancies already found earlier, it is not surprising that 
overall a large spread is also found here (Fig. 114-i . Compared 
to the observations the blackbody SED seems again to fit best. 
The WM-Basic, CMFGEN and CoStar models are marginally 
compatible with the observations each one on one side, while 
the TLUSTY and Kurucz results are definitively far away from 
the observed values. Note, however, that both theoretical quan- 
tities plotted here (including 77) depend also on the ionization 
parameter. 

In Fig. [B] we illustrate theoretical predictions of mid-IR 
softness parameters as a function of T e ff for the case of r\Ne-Ar- 
Note that most of the model atmospheres predict that rj^e-Ar 
becomes insensitive to T e g above a certain value, here of the 
order of T e jf 35000 K (the exact value also depending on 
U). Similar dependencies on the adopted model atmosphere 



Fig. 17. Variation of the radiation softness parameter rjs-Ne ver- 
sus T]A r -Ne, when changing U, T e ff or Z. Arrows codes are the 
same as in Fig. [6] 

and "saturation effects" have also been found for the tr aditional 
optica l softness parameter (cf. lOev et al.l2002t IKennicutt et alJ 

Furthermore model calculations also show that all the mid- 
IR 77 depend quite strongly on the ionization parameter U and 
on metallicity, as shown here for the case of the WM-Basic 
models 7 (see Fig. I17> . From the theoretical point of view, 
and without tailored photoionization modeling including con- 
straints on U and Z, the use of mid-IR 77's appears therefore 
highly compromised. Again, similar difficult ies have also bee n 
found for the o ptical s oftness parameter , cf. Skillman (1989); 
iBresolin et all Jl999h : lOev et al.1 d2002l) . but IKennicutt et alJ 
(2000) and hereafter. 

Last but not least, as each 77 implies 4 line intensities, the 
softness parameters are more sensitive to any observational un- 
certainty (as the attenuation correction or detector calibration), 
possible collisional effects, uncertainties in the atomic data etc. 
In view of all these considerations and the model results pre- 
sented above, applications of softness parameters in the mid-IR 
appear therefore to be of very limited use. 

For completeness we show in Fig. [HO the behavior of the 
traditional optical softness parameter 770-5 8 for all model at- 
mospheres. Again a considerable spread between different at- 
mosphere models and a dependence on U and Z is found. 

7. Discussion 

From the comparisons of the extensive model calculations with 
the observed excitation diagnostics we can draw some rather 
general conclusions on the shape of the ionizing spectra of 

7 Note that the metallicity dependence predicted by the WM-Basic 
models is probably overestimated as already discussed earlier (Sect. 
l5~5l . 

8 Defined bv IVilchez & Paeell 1 19881) as : ([On] 3726,27 A/[0 m] 
4959,5007 A) / ([Sn] 6717,31A/[S m] 9069,9532A) 
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early type stars and the on the nature of the ionizing sources 
of the Galactic H n regions. 

7.1. Implications on stellar SEDs 

Overall the "best fit" SED to the observed excitation diagrams 
(Figs. [3] to |4} was found with blackbody spectra. Does this 
imply that the ionizing fluxes of hot stars are best described 
by the Planck function ? The overall answer is no, but. The 
observations probe (to 1st order) the relative number of ion- 
izing photons with energies above the relevant ionization po- 
tentials, i.e. what was called the "slope" in Sect. 5. Therefore 
if we considered that all observed 3 line ratios are correctly 
reproduced by a blackbody this would imply that the 27.6- 
35.0 and 27.6-41.1 eV slopes (and hence also 35.0-41.1 eV) 
be equal to that of the Planck function (bb) of the same T e g, 

i.e. Q27.6/Q35.0 = Q27.6/Q35.0 and Q27.6/Q41.1 = Q27.fi/Q4u- 
This would therefore represent three "integral" constraints on 
the stellar SED. This is rather strong, but still leaves room for 
the detailed shape of the SED between these energies. Actually 
the agreement between blackbody spectra and [Arm/n] vs. 
[Nem/ii] is better than diagrams involving [Siv/in]. This in- 
dicates that the constraint on the Q27.fi/Q41. 1 is better than at 
intermediate energies. 

However, it is important to remember that these constraints 
on the underlying SED can be deduced only if the observed 
excitation diagram (Figs[3]to[5j is essentially driven by a tem- 
perature sequence (i.e. U = const). Note also the effects of 
uncertainties on the atomic data, as discussed in Sect. 15.71 



7.2. On the origin of the observed excitation 
sequences 

The origin of the observed mid-IR excitation sequences 
and their correlation with galac tocentric dis ta nce ha s 
been discussed recently b y IGiveon et al.l J2002b): 
Martin-Herna ndez et al.l (|25o2b). As both studies present 
somewhat limited arguments a more general discussion is 
appropriate here. The basic question is whether variations of 
the stellar effective temperature or metallicity variations are 
responsible for the observed decrease of excitation toward the 
Galactic Center? 

Give on et al ] (|2002b) have noted from photoionization 
models using WM-Basic spectra that [Ne m/ii] is predicted to 
decrease more rapidly than [Ar in/11] with increasing metallic- 
ity - a finding also confirmed here (cf. Fig.|6j. However, as the 
observed sequence does NOT follow this trend, they conclude 
that the decrease of excitation must be due to a reduction of 
T e ff as opposed to a softening of the stellar SED with increas- 
ing metallicity. Obviously as such this conclusion cannot be 
upheld as the same mid-IR observations, allowing fairly accu- 
rate abundance determin ations, clearly establ ish the existence 
of a metallicity gradient ( IGiveon et alJl2002al) . In fact, the ap- 
parent contradiction between the observed trend with metallic- 
ity and the one predicted with WM-Basic model atmospheres 
indicates quite likely that the stellar SEDs soften too quickly 
with increasing metallicity and/or that the ionization parameter 




2 3 
age (Myr) 

Fig. 18. Average stellar effective temperature and dispersion as 
a function of time for metallicities Z = 0.008 (1/2.5Z , solid 
lines) and 0.04 (2Z , dashed) predicted for an ensemble of sin- 
gle star H 11 regions for a Salpeter IMF with M up = 12OM and 
a minimum Lyman continuum flux log(Qi; m ) > 48.. Note the 
rather small differences in average T e ff and the large intrinsic 
dispersion for each given metallicity. 



in regions at small galactocentric distance must be larger than 
assumed (cf. Sect. l5.5> . 

In c ontrast to the above study, IMartfn-Hernandez et alJ 
d2002bh show a loose correlation between excitation and 
metal abundances (e.g. between [Nem/n] and Ne/H), stress 
the i mportan ce of met allicity effects on the SED (see also 
IMartfn-Hernandez et alJ2004T) . and conclude that at least partly 
the observed decrease of [Ne in/11] must be due to a softening 
of the stellar SEDs with increasing metallicity. 

7.2.1 . Stellar evolution effects 

From what we know, three effects are related to metallicity 
and must all be taken into account. First, higher metallicity is 
known in stellar evolution to lead to a cooler zero age main 
sequence and to an overall shift to cooler temperatures (e.g. 
Sch aller et alJll992l) . Second, blanketing effects in the atmo- 
spheres become stronger with increasin g metallicity and lead 
to softer ionizing spectra (e.g. Sect. 15. 51 Pauldr ach et aljEoOll 



Martfn-Hernandez et al. 2004). Finally, an increased nebular 
abundance leads also to a somewhat lower excitation of the 
gas in the Hn regions (cf. Sect. 15. 5> . The remaining questions 
are then a) which of these effects dominate and b) whether 
taken together they can indeed quantitatively reproduce the en- 
tire range of observed excitation variations in the Galactic H n 
regions. 

Although the content of stellar ion izing sources of the ob- 
jects considered is not known (but see Morisset 120031) we can 
estimate the T e ff variations expected from stellar evolution, 
e.g. by assuming a single ionizing source. We then perform 
Monte Carlo simulations of single star H 11 regions of different 
metallicities assuming that the ionizing stars are drawn from 
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Fig. 19. [Ar m/ii] versus metallicity, measured here by the Ne 
abundance (see text). The line is a linear fit to the observations 
in log-log space. Arrows as in Fig. [6] for WM-Basic models. 
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Fig. 20. [Nem/ii] versus metallicity, as in Fig.^] A very simi- 
lar plot is obtained for [S iv/m]. 



a Salpeter IMF with a given upper mass cut-off M up . In order 
to compute the mean properties of these stars, such as their 
average T e ff, the predicted variation with metallicity and their 
dispersion, the equivalent of a lower mass limit must also be 
specified. This is done by imposing a lower limit on the total 
Lyman continuum photon flux Qn m . Only stars with Q > Qj; m 
at a giv en age are retained fo r this computation. In practice we 
use the iMevnet et alJ {1^94) stellar tracks, we consider metal- 
licities between ~ 1/2 and 2 times solar, as indicated by the 
observed range of Ne/H or Ar/H abundances and we adopt 
log(Qiim) = 48., corresponding a typical lo wer limit for the H n 
regions of lMartfn-Hernandez et al.l J2002al) . Very massive stars 
entering the Wolf-Rayet phase already on the main sequence 
are also excluded. 

The resulting average and spread of T e ff as a function of age 
is shown in Fig.^Jfor Z = 0.04 and 0.008 respectively. This 




figure shows the following: first, the reduction of the average 
stellar temperature due to a metallicity increase by a factor 4 
is rather modest, of the order of ~ 3-4 kK. Second, for a given 
metallicity the predicted dispersion in T e ff is larger; typically of 
the order of ~ 3-9 kK for reasonable ages. Although the abso- 
lute values of these T e ff depend for obvious reasons on the exact 
choice of Qy lm , the T e ff differences and dispersion depend little 
on this value. From all the models considered above, the de- 
crease of the mean T e ff due to stellar evolution effects appears 
to be too small to explain the full range of excitations. 

On the other hand the fairly large T e ff dispersion at a given 
Z will induce quite important variations in the excitation. This 
effect probably dominates the observed excitation scatter (de- 
fined at the end of Sect. [5}, as also suggested by the observa- 
tions of a large spread of the excitation diagnostics Figs.ll9land 
a given metallicity measured here by Ne/Ne Q . 

Indeed the observed spread of [Arm/ii] and [Nem/ii] at a 
given Z is somewhat larger than or similar to the decrease of the 
mean excitation with increasing Z. Taken together these find- 
ings imply that while undeniably metallicity effects on stellar 
evolution and nebular abundances must be present, statistical 
fluctuations of the effective stellar temperature due to the IMF 
are likely the dominant source of scatter for the observed mid- 
IR excitation sequence of Galactic H n regions, while the exci- 
tation sequences must be predominantly driven by other effects 
which we will discuss now. 

7.2.2. Stellar atmospheres and nebular effects 

As apparent from our modeling (see Figs. I19land 120b the ef- 
fects of metallicity on the shape of the stellar ionizing spectra 
strongly alter the predicted excitation. The magnitude of the 
predicted effect is found to be comparable to the observed vari- 
ation. Both these findings and the above results concerning stel- 
lar evolution effects indicate the Z dependence of the ionizing 
spectra is the main driver for the correlation of the excitation 
with galactocentric distance. 

As this result depends on a specific set of model atmo- 
spheres (WM-Basic) a few words of caution are, however, nec- 
essary here. First, we note that the effects of T e tf and Z are 
not exactly the same for [Arin/n],[Ne m/ii], and [Siv/m]. The 
predicted excitation variation (shown here for a change of T e ff 
from 35 to 40 kK) tends to be somewhat larger (smaller) for 
[Nem/n] and [S iv/m] ([Arm/n]) than the observed variations. 
Second, it must be remembered that the WM-Basic atmosphere 
models employed here could predict too strong a softening with 
increasing Z as suspected from Fig.|20|and already discussed in 
Sect. 15.51 Despite these imperfections there is little doubt that 
the above result remains valid. 

Finally we may also comment on excitation changes related 
to the ionization parameter. Again, as for T e ff, the above Monte 



9 Ne abundanc es were obtained using reddening and T e i ec corrected 
abundances from lGiveon et alJ l2002ah . Ne is preferred to Ar or S be- 
cause Ar abundance is not reliable for high T cff , Ar 3+ being present 
but unobserved, the same applying at low T eff for S + . The solar abun- 
danc e for Ne is deter mined using the abundance gradient obtained by 
iGiveon et al]j2002al) at 8.5 kpc. 
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Carlo simulations show small differences between the average 
ionizing photon flux Qn.f, with metallicity, but a considerable 
dispersion for each Z. Combined with the observational fact of 
fairly similar nebular densities in our objects this could be a 
justification for a constant ionization parameter, at least on av- 
erage. Random variations of U are, however, expected to con- 
tribute to the excitation scatter at a given metallicity. 

In conclusion we see little doubt that the observed excita- 
tion sequence of Galactic H n regions is shaped by the joint ef- 
fects of metallicity on stellar evolution, atmospheric line blan- 
keting, and cooling of the ISM. From our investigations it 
seems, however, that metallicity effects on ionizing stellar flux 
is the dominant effect causing the excitation gradients while 
statistical fluctuations of T e ff and U are likely the dominant 
source of scatter in the observed excitations. A more detailed 
study of possible sys tematic TW and ionization parameter gra- 
dients is presented in Morisset (2003): no clear gradient of T e g 
nor U versus the galactocentric distance are found by this au- 
thor, but some trends of increase (decrease) of T e ff (U) with 
the metallicity are observed. Importance of taking into account 
the effect of the metallicity on the stellar spectral shape is ad- 
dressed. 



8. Summary and Conclusion 

We have presented results from extensive photoionization 
model grids for single star H n regions using a variety of recent 
state-of-the-art stellar atmosphere models such as CMFGEN, 
WM-Basic, TLUSTY, CoStar, and Kurucz models. In fact, even 
among the two recent non-LTE line blanketed codes including 
stellar winds (WM-Basic and CMFGEN) the predicted ionizing 
spectra differ by amounts leading to observable differences in 
nebular spectra 10 

The main aim of this investigation was to confront these 
model predictions to recent catalogs of ISO mid-IR observa- 
tions of Galactic H n regions, which present rich spectra prob- 
ing the ionizing spectrum between ~ 24 to 41 eV thanks to 
the measurements of [Ariii/n], [Nem/n], and [Siv/m] line ra- 
tios. Particular care has been paid to examining in detail the 
dependences of the nebular properties on the numerous nebu- 
lar parameters (ionization parameter U, abundances, dust etc.) 
which are generally unconstrained for the objects considered 
here. 

Most excitation diagnostics are found to be fairly degen- 
erate, but not completely so, with respect to increases of T e g, 
U, a change from dwarf to supergiant spectra, a decrease of 
the nebular metallicity (Sects. l5~2l and l5~3l . and the presence 
of dust in the Hn region (Sect. f5~6t . Each of these parameters 
increases the overall excitation of the gas, and in absence of 
constraints on them, a derivation of such a parameter, e.g. an 
estimate of the stellar T e ff of the ionizing source, is intrinsi- 
cally uncertain. In consequence, while for sets of objects with 
similar gas properties statistical inferences are probably mean- 



10 Presumably these differences are due to the use of different meth- 
ods to treat line blanketing (opacity sampling method versus super- 
level approach in the comoving frame) and different atomic data. 



ingful, such estimates for individual objects must be taken with 
care. 

Provided the ionization parameter is fairly constant on av- 
erage and the atomic data is correct (but cf. below) the compar- 
isons between the photoionization model predictions and the 
observations allow us to conclude the following concerning the 
different stellar atmosphere models (Sect. l5"Tt : 

- Both recent non-LTE codes including line blanketing and 
stellar winds (WM-Basic and CMFGEN) show a reason- 
able agreement with the observations. Given their different 
behavior in the three excitation diagnostics, depending on 
luminosity class, and other remaining uncertainties, it ap- 
pears that none of the models can be preferred on this basis. 

- The plane parallel hydrostatic codes (Kurucz, TLUSTY) 
predict spectra which are too soft, especially over the en- 
ergy range between 27.6, 35.0, and 41.1 eV and above. 
Although a good agreement is found for UV to opti- 
cal spectra predicted by the hydrostatic TLUSTY code 
and the photosphe re-wind code CMFGEN jHillierl 12003: 
Bour et et alj E003) important differences are found in the 
EUV range probed by the present observations and pho- 
toionization models. Apparently the full non-LTE treatment 
of numerous elements accounted for by TLUSTY is insuffi- 
cient to accurately predict the ionizing spectra at these en- 
ergies, and the inclusion of stellar winds is imperative. 

- We confirm th e finding of earlier investigations (e.g. 
lOev et alll2000l) showing that the CoStar models overpre- 
dict somewhat the ionizing flux at high energies. 

- Interestingly blackbodies reproduce best the observed ex- 
citation diagrams, which indicates that the ionizing spectra 
of our observed object should have relative ionizing pho- 
ton flux productions Qe at energies 27.6, 35.0 and 41.1 eV 
close to that of blackbody spectra. Although this integral 
constraint on the SED remains approximate, it should still 
be useful to guide future improvements in atmosphere mod- 
eling. 

- Finally, the softening of the ionizing spectra with increas- 
ing metallicity predicted by the WM-Basic models is found 
to be too strong. As already apparent from observed corre- 
lations between excitation diagnostics probing various en- 
ergies, the observed softening of the radiation field (in part 
due to metallicity) affec ts fairly equally the ran ge between 
~ 27 and 41 eV ( M artin-Hernandez et al J2002al) in contrast 
to the atmosphere model predictions, which soften most at 
the highest energies. 

These conclusions are found to be fairly robust to effects 
such as changes of U, nebular and stellar metallicity changes, 
and the inclusion of dust. We suggest that the main uncertainty 
which could alter the above conclusions is the poorly known 
atomic data for Ar ++ (especially dielect ronic recom b inatio n 
coefficients) as also pointed out by Stasihsk a et alJ (^002). 
Reliable computations for such data are strongly needed. From 
the perspective of atmosphere codes probably the most impor- 
tant step toward improving the reliability of ionizing fluxes re- 
sides in a quantitative exploration of the influence of X-rays on 
the emergent spectra at lower energy. 
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The potential of mid-IR line ratios or "softness parame- 
ters", defined in analogy to the well known 77 parameter for 
optical emission lines, has been explored (Sect. loTt . The fol- 
lowing main results have been obtained: 

- Given the non-negligible differences between the various 
atmosphere models it is not surprising that individual line 
ratios (e.g. [Arm/ii], [Nem/ii]) show quite different depen- 
dences on T e ff. We find that [Arm/11] depends little on the 
ionization parameter as the ionization of Ar + is closely cou- 
pled to that of He. This suggests that [Arm] 8.98/-/m/[Arii] 
6.98/vm should in principle be a fairly robust temperature 
indicator, provided the atmosphere models are sufficiently 
accurate up to ~ 24-27 eV and the atomic data is reliable 
(cf. above). 

In comparison to He i/H indi cators (e.g. H e i 667 8A/H/3 or 
Hei 2.06yum/Bry in iKennicutt et aljEoOdt iLumsden et alJ 
2003, and references therein), [Arm/11] does not show a 
saturation effect but remain sensitive to T ef j up to the high- 
est temperature examined here (T ef f ~ 50 kK). Due to the 
large uncertainties of die lectronic recombination coefficient 
of Ar ++ , Morisset ( 2003) prefer to use [S iv/m] and [Ne m/11] 
to determine T e ff and U simultaneously, for the H 11 regions 
used in this work. 

- Both empirically and theoretically the mid-IR softness pa- 
rameters which can be constructed from [Arm/ii], [S iv/m], 
and [Nem/ii] are found to provide little if any informa- 
tion on stellar temperatures i f not used to determine U at 
the same time (Morisset 2003). Observationally little / no 
correlation is found between the 77's and excitation ratios. 
Furthermore, our photoionization models show a consider- 
able dependence of 77 on the ionization parameter. We there- 
fore conclude that mid-IR 77's appear to be of limited diag- 
nostic power. 

Finally we have examined which parameter(s) is (are) 
chiefly responsible for the observed mid-IR excitation se- 
quences. Combining the results from our extensive photoion- 
ization model grids with Monte Carlo simulations of ensem- 
bles of single star H 11 regions of different metallicity and age 
we conclude the following (Sect. 0. While metallicity effects 
on stellar evolution, atmospheres and the nebulae all have an 
undeniable influence, they are probably of minor importance 
compared to the fairly large dispersion of T e ff expected at each 
metallicity from a simple statistical sampling of the IMF. The 
T e ff scatter plus additional scatter in the ionization parameter 
are probably the dominant driver for the observed mid-IR exci- 
tation scatter of Galactic H 11 regions, while the effect of metal- 
licity on the shape of the ionizing spectra is partially respon- 
sible of the global excitation sequence, the proportion of this 
effect being str ongly dependen t of the reliability of the atmo- 
sphere models jMorisse 120031) . 
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